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Abstract 

We compute several coefficients needed for 0(a) improvement of currents in perturbation the- 
ory, using the Brodsky-Lepage-Mackenzie prescription for choosing an optimal scale q* . We then 
compare the results to non-perturbative calculations. Normalization factors of the vector and axial 
vector currents show good agreement, especially when allowing for small two-loop effects. On the 
other hand, there are large discrepancies in the coefficients of 0(a) improvement terms. We suspect 
that they arise primarily from power corrections inherent in the non-perturbative methods. 
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I. INTRODUCTION 



During the past few years the Symanzik effective field theory has been an important focus 
of research in lattice gauge theory. Symanzik's idea is to describe cutoff effects in lattice 
field theory by a continuum effective field theory [|lj]. One writes [|I|, 0] 

Aat = ^Sym, (1) 

where the symbol = means that the lattice and Symanzik field theories have the same on- 
shell matrix elements. For lattice QCD with Wilson fermions the Symanzik local effective 
Lagrangian (LE£) is given by [|, |J 

£ Sym = 5-2 tr[F^F^] - q (p + m) q + aK a . F qia^F^q + 0(a 2 ), (2) 

where g 2 is a renormalized coupling, m is a renormalized quark mass, and aK a .p is a short- 
distance coefficient. The effective field theory is useful when the scale of QCD in lattice units 
is small, Aa <C 1, and, as used in this paper, when ma <C 1 also. With the description in 
hand, the lattice field theory can be adjusted so that it approaches its continuum limit more 
quickly. The effective theory shows that if K a .p is reduced for any given on-shell matrix 
element, then the 0(a) term in Eq. (0) makes commensurately smaller contributions to all 
other on-shell matrix elements. This application of the Symanzik effective field theory is 
called the Symanzik improvement program ||. 

A similar correspondence is set up for the vector and axial vector currents (see below), 
introducing further short- distance coefficients. In the last several years methods have been 
devised to study all of them non-perturbatively |6|, [7|, |], H |10fl - The 0(a) discretization ef- 
fects violate chiral symmetry, so the key idea is to ensure that violations of chiral symmetry 
are at least 0(a 2 ). On the other hand, because of asymptotic freedom and the success of 



perturbative QCD, even at GeV energies [11], one expects perturbation theory to yield ac- 



curate estimates of the short- distance coefficients. In this paper, we compare a perturbative 
calculation of the currents' short-distant coefficients to the non-perturbative results. 

There are two issues that should be kept in mind when making such a comparison. First, 
the non-perturbative technique suffers from power corrections. Asymptotically, as Aa — > 
these are formally smaller than any error made from truncating the perturbation series. In 
practice, however, these effects can be significant. 

Second, no two-loop results are available for the improvement coefficients considered here. 
Tests of perturbation theory are, therefore, not unambiguous, because different choices for 
the expansion parameter g 2 yield quantitatively different results. The bare coupling g 2 (for 
the Wilson gauge action) is an especially bad expansion parameter [ F2~ . The obvious remedy 



is to rearrange the perturbative series, eliminating g 2 in favor of a renormalized (running) 
coupling, evaluated at a scale characteristic of the problem at hand. One is then faced, 
however, with many choices of renormalization scheme, and the question of how to determine 
the "characteristic scale." In this paper we choose the Brodsky-Lepage-Mackenzie (BLM) 



prescription [13|, p4| . Once this choice is made, little subjectivity remains, so one can ask 
quantitatively whether one-loop BLM perturbation theory agrees with the non-perturbative 
method. 

In the BLM method, the characteristic scale is computed from Feynman diagrams. The 
new information presented in this paper consists of the calculations needed to determine the 
BLM scales of the normalization and improvement coefficients of the vector and axial vector 



2 



currents for Wilson fermions with Sheikholeslami-Wohlert action. These calculations are 
a by-product of our recent work on the normalization and improvement of lattice currents 
with heavy quarks fI5 |. Details of the calculational method may be found there. 



This paper is organized as follows. In Sec. [Tl]we define the lattice currents and review their 
description in the Symanzik effective field theory. Section [TT1] recalls the BLM prescription, 
focusing on points that are sometimes overlooked. Our new results for the BLM scales are 
given in Sec. [TV]. This paves the way for a systematic comparison with non-perturbative 
calculations of the same quantities in Sec. |V|. Section |VJ contains a few concluding remarks. 



II. LATTICE CURRENTS 

In this section we review the description of lattice currents with the Symanzik effective 
field theory. For quarks we take the Sheikholeslami-Wohlert Lagrangian [|]], which has 
an improvement coupling csw- At the tree level K^} F = |(1 — csw), so the improvement 
condition K a .p = requires csw = 1 + 0(g 2 ). For one-loop calculations, it is sufficient to 
specify csw at the tree level. For the non-perturbative calculations cited below, csw — 1 is 
determined non-perturbatively by the methods of Ref . . 

We denote the lattice fermion field with if>. The lattice vector and axial vector currents 
take the form 

V& = fa-fiff - ac v d vl ^a^, (3) 
A L = iri-'flfti + ac A d^ t ipij 5 ip. (4) 

The improvement couplings cy and ca should be chosen to reduce lattice artifacts, as dis- 
cussed below. 1 In Symanzik's theory of cutoff effects, the lattice currents are described by 
operators in a continuum effective field theory [|l], |J [5], [/]] 

V& = Zy l qi^q - aK v d v qa^q + • • • , (5) 
Af at = Zfqi^q + aK A d»qi l5 q + ■ • ■ , (6) 

where, as in Eq. (^), q is a continuum fermion field whose dynamics is defined by jCqcd- The 
ellipsis indicates operators of dimension five and higher. Further dimension-four operators 
are omitted from Eqs. ([5]) and (|6|), because they are linear combinations of those listed and 
others that vanish by the equations of motion. The short-distance coefficients in the effective 
Lagrangian — Zj and Kj (J = V, A) — are functions of g 2 and ma, and the improvement 
couplings csw and c j. 

Symanzik improvement is achieved by adjusting cj so that Kj = 0. Then ZyV^ t and 
ZaA^ have the same matrix elements as qi^q and qi'j^'j^q, apart from lattice artifacts of 
order a 2 . For light quarks one may expand Zj in ma, 

Zj = Zj{l + mabj), (7) 



1 The lattice currents in Eqs. (Q) and (TJ) are useful for light quarks. For heavy quarks the "small" improve- 
ment terms become large, introducing unnecessary violations of heavy-quark symmetry. Better currents 
for heavy quarks are given in Refs. |]l5[ |l6| |. 
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and identify Kj with the zeroth order of a small ma expansion. At the tree level the 
coefficients of the normalization factor are Zj' = 1, = 1. In addition, the coefficient of 
the lattice artifact is 

Kf = cM (8) 

The improvement condition Kj = says that one should set Cj = 0. Consequently, one-loop 
calculations are based solely on the first terms in Eqs. (0) and (f|). 



III. BRODSKY-LEPAGE-MACKENZIE PRESCRIPTION 

In this section we review the BLM prescription, following the argumentation from 



Ref. [14]. This material should be familiar, but some of the literature on non-perturbative 
improvement blurs the difference between BLM perturbation theory and other topics, such 
as "tadpole improvement" and mean-field estimates of the renormalized coupling, which are 



also discussed in Ref. [14 



The problem is to find a reasonably accurate one-loop estimate of a quantity £, here Zj 
or Kj. In these cases, one gluon with momentum k and propagator D(k) appears. The 
contribution from the Feynman diagrams can be written 

9 2 R C [1] (p)=9oj j0jiD{k)f{k,p) + ..., (9) 

where p denotes fc-independent parameters, such as external momenta. The ellipsis indicates 
higher-order terms that we would like to absorb into the renormalized coupling g\. An 
important class of higher-order terms consists of the renormalization parts that dress the 
exchanged gluon. In the Fourier transform of the heavy-quark potential, for example, they 
turn g^D(k) into g v {k)D{k), where the potential V(q) = —CFg v (q)/q 2 . Thus, 

9%C [1] (P) = J 0^g 2 v (k)D(k)f(k,p) + ■■■ (10) 

sums the renormalization parts. Other ways of dressing the gluon would lead to other 
physical running couplings, but they all are the same at order /3 g 4 where j3 = 11— |n/ 
is the one-loop coefficient of the (3 function for rif light quarks. 
If there is a characteristic scale q*, one can approximate 

fl 2/i.\ = 9y((l*) (u) 

gv[ ' l + ((3 /W^)g v (q*)Hk/q*y 1 ) 

= gW) + ^- 2 9y{q*)Hq*/k) 2 + ---. (12) 

The aim is to choose q* so that higher-order terms are small, particularly those of order 
/3 <7y, which could be enhanced by a foolish choice of q*. Inserting Eq. (|T2| ) into Eq. ( [TDD 
and setting the coefficient of Po9v t° zero yields 

*£[i] 

lng * a= 2CW' ( 13 ) 
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where a is a reference short-distance scale (namely, the lattice spacing), and 



xll]{p) = I ^M^m^M- ( 14 ) 

Thus, the BLM prescription is to set g 2 R = g v (q*) in the one- loop approximation. 

If one prefers a different renormalized coupling, one must change the scale in the appro- 
priate way. The coupling in scheme "S"' is related to the V scheme by 

+ P ° b "\ bs +0(g% (15) 



where and bg are constants independent of n/. The BLM scale q* s for this scheme is 
given by 

lng* =\nq*-±bf. (16) 
For example, for the modified minimal subtraction (MS) scheme, b^L = —8 and b^ = 5/3, 

1V1 o 1V1 o 

^Ms = e_5 ^ 6 <?* = 0.435g*. With Eq. (0) one recovers the summary statement of Ref. Ul3| , 
namely to absorb into the rif dependence of the two-loop term, which enters through (3q. 

The BLM prescription has several features that make it a natural choice in matching 
calculations, such as those considered in this paper. The effective field theory framework 
suggests using a renormalized coupling, in particular one that has a (quasi-)physical defi- 
nition in both the underlying theory (here lattice gauge theory) and in the effective theory 
(here the Symanzik effective field theory). For quantitative purposes it is more interesting 
to note that 

1 1 b {0) 

+ T^ + Otf), (17) 



9 2 s(q*s) 9v(q*) 16vr 2 

so the numerical difference in the BLM expansion parameters is small, as long as g 2 b^ /167T 2 
is small. 



IV. PERTURB ATIVE RESULTS 



In Ref. |15|] we found for gauge group SU(3) and csw = 1 

Z [ y ] = -0.129423(6), (18) 
Zj? = -0.116450(5), (19) 



in excellent agreement with previous work [|17|, |18| . (Reference |18| gives precise results as a 
polynomial in csw-) We also found (with Cj 1 = 0) 

b v ] = 0.153239(14), (20) 

b [ l ] = 0.152189(14), (21) 

K [ y ] = c v ] +0.016332(7), (22) 

K [ } ] = cJJ 1 +0.0075741(15), (23) 
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which agree perfectly with Ref. [19|. Solving the improvement condition Kj = at this 
order gives 

C W = -0.016332(7), (24) 

cJJ 1 = -0.0075741(15). (25) 

We also directly obtained 

&M _ 6 W = 0.0010444(16), (26) 

which is more accurate than the difference of the two numbers quoted above. In taking the 

difference, large contributions from the self energy cancel, but, even so, the near equality 

of by and is a bit astonishing. The mass dependence of Zj shows that fry — is not 
so small for the Wilson action ||15|| . 

In our method for computing the improvement coefficients it is easy to weight the inte- 
grands with \n(ka) 2 and, thus, obtain the BLM scales. We find 

*Z [ fi = -0.270691(19), (27) 

*Z [ j ] = -0.243086(09), (28) 

*b v ] = 0.321556(35), (29) 

*b [ A ] = 0.318108(21), (30) 

*&W-*&5J ] = 0.0034247(51), (31) 

* C W = -0.0222383(15), (32) 

*c[} ] = -0.0147825(62), (33) 

and hence 



q* Zv a = 2.846, (34) 

q* ZA a = 2.840, (35) 

q* ZA/Zv a = 2.898, (36) 

q* bv a = 2.855, (37) 

q* bA a = 2.844, (38) 

ql v _ bA a = 5.153, (39) 

q* cy a = 1.975, (40) 

q* CA a = 2.653. (41) 

The scales are in the expected range. The higher scale for by — bA means simply that 
the difference between these renormalization constants arises from very short distances. 
These numerical results are new; they have been obtained from two independent computer 
programs. As a further check, we have reproduced the values of q* Zv a and q* ZA a for the 
Wilson action (csw — 0), given in Ref. |20| . 



The dominant contributor to the "large" one-loop normalization constants, Eqs. 

is the tadpole diagram (in Feynman gauge) of the self energy. One might expect 
perturbation theory to work better for quantities in which the effects of tadpole diagrams 
largely cancel (albeit in a gauge-invariant way). For example, Za/Zv and by — bA are tadpole 
free and have smaller one-loop coefficients. 
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Another way to remove the tadpoles is to write 



Zj = u Zj, (42) 
bj = h/uo, (43) 

where uq is any convenient tadpole-dominated quantity Then one can take uq from a 
non-perturbative Monte Carlo calculation and use perturbation theory for Zj and bj. The 
corresponding one-loop coefficients are 

Zf = Zf - «W (44) 

Sg" = tS" +0 W. (45) 

Similarly, to get the BLM scale 

*Z l j ] = *Z [ j ] - *v$\ (46) 

*b [ j ] = *bf + *v$\ (47) 

where *Uq is the BLM numerator [cf. Eq. fll3|) 1 for Mq. Below we take Uq to be the average 
value of the plaquette, with hJ 1 = -1/12 = -0.083 and *uf = -0.204049(1). A glance 
at Eqs. (|44])-(|47D shows immediately that tadpole improvement reduces the one-loop coeffi- 
cients. With tadpole improvement the BLM scales become 

q\ y a = 2.061, (48) 

q\ A a = 1.803, (49) 

qi bv a = 2.317, (50) 

q\ A a = 2.289. (51) 

The scales are lower than without tadpole improvement, but still ultraviolet. 

It is perhaps worthwhile emphasizing the difference between tadpole improvement and 
the BLM prescription. The aim of tadpole improvement is to re-sum large contributions 
appearing at order g 2 and higher, replacing the sum with a non-perturbative estimate (uq, 
for example). The aim of the BLM prescription is to re-sum potentially large renormalization 
parts into the renormalized coupling. Although the aims are similar, they are not identical. 
They are not mutually exclusive, and neither is a substitute for the other. 



V. COMPARISON TO NON-PERTURBATIVE CALCULATIONS 



With the BLM scales in hand we can compare the prediction of one-loop BLM-improved 
perturbation theory with non-perturbative determinations of the improvement coefficients. 
We shall make the comparison in two ways. First we compare the numerical values directly, 
at two values of the bare coupling. Here there are two methods in the literature, one 

and another 



based on finite-size techniques and the Schrodinger functional [El 



23 



based on large volumes with hadronic matrix elements ||24|| . The difference between these 
two illustrates how large power corrections to the improvement coefficients are. We also 
compare our results graphically, as a function of coupling, to Pade approximants given in 
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TABLE I: Comparison of perturbative and non-perturbative determinations of the improvement 
coefficients at (3 = 6.2. 



(3 = 6.2 


otv ( o*) 


BLM 


Refs. |2T], |23| 


Ref. fl2l 


Zv 


0.f468 


0.7612 


7922(9) 


7874(4) 


Za 


0.1469 


0.7850 


0.807(8) 


0.818(5) 


Z A Z\f 


0.1461 


1.0238 


1.019(8) 


1.039(5) 


bv 


0.1467 


1.2824 


1.41(2) 


1.42(1) 


I) A 


0.1468 


1.2808 




1.32(5) 


by ~ bA 


0.1257 


0.001649 




0.11(5) 


-c v 


0.1638 


0.03361 


0.21(7) 


0.09(2) 


-CA 


0.1498 


0.01426 


0.038(4) 


0.032(7) 


u Z v 


0.1616 


0.8022 


0.7922(9) 


0.7874(4) 


uqZ a 


0.1686 


0.8230 


0.807(8) 


0.818(5) 


bv/uo 


0.1559 


1.2846 


1.41(2) 


1.42(1) 


b A /u 


0.1565 


1.2828 




1.32(5) 



TABLE II: Comparison of perturbative and non-perturbative determinations of the improvement 
coefficients at f3 = 6.0. 



P = 6.0 


a v (q*) 


BLM 


Refs. H 


Ref. H 


Z v 


0.1602 


0.7394 


0.7809(6) 


0.770(1) 


Z A 


0.1603 


0.7654 


0.791(9) 


0.807(8) 


Za/Zv 


0.1593 


1.0260 


1.012(9) 


1.048(8) 


bv 


0.1601 


1.3082 


1.54(2) 


1.52(1) 


b A 


0.1603 


1.3065 




1.28(5) 


bv ~ bA 


0.1352 


0.001774 




0.24(5) 


-c v 


0.1808 


0.03711 


0.32(7) 


0.107(17) 


-CA 


0.1638 


0.01559 


0.083(5) 


0.037(9) 


u Z v 


0.1782 


0.7872 


0.7809(6) 


0.770(1) 


u Z A 


0.1868 


0.8095 


0.791(9) 


0.807(8) 


bv/uo 


0.1712 


1.3105 


1.54(2) 


1.52(1) 


b A /u Q 


0.1719 


1.3087 




1.28(5) 



Refs. ||19| , Em . These graphs are helpful for seeing whether discrepancies in the one-loop 
and non-perturbative estimates arise from two-loop or power corrections. 
We obtain «y(g*) as follows. First we compute 



«lxl 



(52) 



where (□) is the ensemble average of the plaquette. Then we follow Ref. and take a>v 
to be 

«"("02/a) = 1 + v ,^, OM (53) 
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FIG. 1: Zy and Za vs. g$. Heavy lines show the non-perturbative results, Eqs. (|54f) and (|55|), and 



shading possible corrections of order ±(Aa) 2 . Circles show BLM perturbation theory, with thin 
and dashed lines to indicate a two-loop term ia^ or ±2a v . (a) and (c) no tadpole improvement, 



^BLM 



1 + 9v(<l*Zj} Z J ] i ( b ) and ( d ) with tad P° le improvement, Z^ LM = u [l + g v {q\ )Z 



r[l]l 



which agrees with the standard definition of ay with an accuracy of order o? s . The scale 
3.402/a is the BLM scale for (□). We then run from 3.402/a to q* with the two-loop 
evolution equation. Of course, once two-loop perturbation theory is available, one would 
have to extend the accuracy of Eq. (|53"D and of the evolution. 

Table HI gives results from our perturbative calculation with non-perturbative results from 
the Alpha Collaboration |2l], |22|, and from Bhattacharya et al. ||24|| , at (3 = 6.2. Table |T| 
gives the same comparison at (3 = 6.0. Above (below) the horizontal line, we have applied 
the BLM prescription without (with) tadpole improvement. The error bars on the entries 
from Refs. |21~| , |22| , p3 , p4]l are statistical, and compiled in Ref. ||24j| . 

For the normalization factors Zy and Z4, BLM perturbation theory and the non- 
perturbative methods agree well, within 3-4%. The difference between the two non- 
perturbative values of Zy exceeds the reported errors, but is easily explained by power cor- 
rection of order (Aa) 2 . For the tadpole-free ratio Z^/Zy and for the tadpole-improved quan- 
tities uqZj, BLM perturbation theory lies very close to the non-perturbative range. These 
impressions are strengthened by Fig. [I], which shows Zy and Za as functions of g$. Circles 
show BLM perturbation theory, and the thin solid (dashed) lines indicate how two-loop 
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FIG. 2: by vs. <7g. (a) no tadpole improvement, b v LM = 1 + g v {ql )by , (b) with tadpole improve- 
ment, 6y LM = [1 +gy(q l T )b^y]/uo. Light grey (light blue) shading indicates power corrections to 
by of order ±Aa; darker grey (pink) shading power corrections to by — 1 of order ±a/L. 
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contributions of ±a y (±2a v ) could modify the result. We show the result with and with- 
out tadpole improvement in Figs. |l](b,d) and (a,c), respectively. For the non-perturbative 
method, a heavy (blue) line shows the Pade approximants 



= 1 - 0.7663^ + 0.048% 4 
V 1 - 0.6369^ ' 1 ' 

l- 0.8496g 2 + 0.0610^ 
A l- 0.7332^ 2 ' 1 ' 

which deviate from the underlying calculations negligibly for g 2 < 1. The shaded bands 
behind the Pade curves show a power-correction of ±(Aa) 2 , with A ~ 500 GeV. The finite- 
volume result also suffers from power corrections of order (a/L) 2 . They are estimated to be 
small by comparing calculations on lattices with a/L = 1/8 and 1/12 Also, they are 
parametrically smaller, because Ref. [22] holds LA ~ 2 for all g 2 . 

Next let us turn to the 0(ma) corrections to the normalization factors, by and b^- There 
is only one calculation of b^ p4j, so let us concentrate first on by. The two non-perturbative 



results for by agree perfectly with each other (see the Tables), but they deviate significantly 
from one-loop BLM perturbation theory. Some insight can be gleaned from Fig. 0, which 
shows by as a function of g 2 . The non-perturbative method is represented with the Pade 



approximant |05 



1 - O.76130Q 2 + 0.00124 4 - 0.1136^ 6 
W ~ 1 - O.91450 O 2 ' (56) 

with light (blue) shading for a power correction ±Aa. In finite volume there is also a power 
correction to by of order a/L; by construction it applies to by — 1 p2[ , but now L with a 
varies such that a/L = 1/8 for all g 2 . We model this effect as {by — 1)(1 ± |), shown in the 
darker (pink) shading in Fig. |2|. Judging from its size and shape, the deviation seen in Fig. |2| 
looks less like a two-loop effect than a combination of power corrections of order a/L and Aa. 
(Similar conclusions are reached in Ref. |24j].) There is almost no difference whether one 
applies tadpole improvement to by or not, once the BLM prescription is applied. These two 
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approximations truncate higher orders of the perturbation series differently, substantiating 
the idea that the discrepancy is a power correction. 

The non-perturbative calculation of 6 a agrees with one-loop BLM perturbation theory. 
Note, however, that Ref. [^4]] obtains by and by — bA directly, and then b& = by — {by — 6,4). 
The agreement between BLM perturbation theory and the non-perturbative results for by 
and by — bA is not good, so the agreement for 6^ may be an accident. Since the coefficient 
by — b^ in Eq. ( [ZED * s remarkably small, the two-loop contribution could be as large as the 
one-loop term. Furthermore, inspection of Fig. 14 in Ref. [El] suggests that a fit to the three 



smallest masses would yield a smaller value of by — bA- We consider the comparison of bA and 
by — bA to be unsettled pending a two-loop calculation and a more robust non-perturbative 
calculation. 

In any case, the mild disagreement on by and by — bA is not of much practical importance. 
For the sake of argument, suppose ma < 0.1, which holds for the light quarks for which the 
currents were designed. Then power corrections in bj, at fixed a, lead to an uncertainty in a 
decay constant or a form factor of only a few per cent. After a continuum limit extrapolation, 
these uncertainties will not be important. 

Now let us turn to the coefficients cj of the improvement terms in Eq. (0) and @. At 
the tabulated values of f3, the non-perturbative and BLM calculations of ca do not agree at 
all. At /3 — 6.0 (Table H) the two non-perturbative calculations also do not agree with each 
other. Figure |3](a) shows function of using the Pade approximant pl| 



,1 - 0.748^ , . 

C * = -°-° 0756s °T^77t (57) 

to represent the non-perturbative calculations. The disagreement between BLM perturba- 
tion theory and Eq. fl57|) sets in for g% > 0.9. There are two reasons to suspect that the 
discrepancy stems from a power correction of order Aa to the results of Ref. [pl| . First, 



Fig. y(a) shows that it has the shape and size of such a power correction. Second, the ex- 
tracted value of ca depends on the lattice derivative used to define the current p5[] . Note f24 



that errors in ca propagate to cy, because in the Ward identities ca is multiplied by large 
hadronic matrix elements such as am 2 K /m s ~ a x 2.5 GeV. This enhancement also explains 
why Eq. (|57D leads to worse scaling in /„- p5| . Figure 0(a) also includes the non-perturbative 
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TABLE III: Expansion parameters for perturbation theory. 



P "0 aixi <*v(q*z v ) a Ms(gk-) 

6.0 0.0796 0.1340 0.1245 0.1602 0.1784 

6.2 0.0770 0.1255 0.1166 0.1468 0.1619 

6.4 0.0746 0.1183 0.1101 0.1362 0.1491 

7.0 0.0682 0.1016 0.0951 0.1134 0.1222 

9.0 0.0531 0.0702 0.0667 0.0748 0.0786 



results of Refs. [^3, F25|. The difference between those points and BLM perturbation theory 
could be a modest two-loop effect or a small power correction. 

For cy, the two non-perturbative results agree neither with each other, nor with BLM 
perturbation theory. The Alpha Collaboration has only a preliminary calculation p3| . We 



have taken the liberty of extracting results from Fig. 3 of Ref. |23| and fitting them to a 
Pade formula. The leading behavior is fixed to Eq. (p4|), and we obtain 



c v = -0.0163340 2 1 °- 257 ^ . (58) 
v y ° 1 - 0.963^ V ; 

Figure |3|(b) plots Eq. fl5"8"D, the underlying points P^j , the non-perturbative results from 
hadronic correlation functions [24"], and BLM perturbation theory. As usual we show possible 



power corrections to Eq. (^) of order ±Aa, as well as the size of typical two-loop effects. 
At small #q, there is good agreement with (BLM) perturbation theory, but once g$ > 0.9, 
there is a sharp turnover. It is probably a power correction, possibly exacerbated by power 
corrections to ca as modeled by Eq. ([57]). With hadronic correlation functions p4| the non- 



perturbative value of cy is half or a third as large. It is not clear at present whether the 
discrepancy between Ref. |24[ and BLM perturbation theory is a power correction to the 



former or a sizable two-loop correction to the latter. 

We should also mention that BLM perturbation theory works better than several forms 



of mean- field perturbation theory (let alone bare perturbation theory). In Table |T| we list 
several choices for a s : 

a = gl/Air, (59) 

a = a /u 0} (60) 

as well as «i x i [Eq. ([52])] and aMs(^Ms) With only one-loop expansions avail- 

able, the mean- field choices «o an d «ixi give smaller corrections, and one- loop perturbation 
theory falls short even when power corrections are negligible. The consistency of BLM-V 
perturbation theory for Z v , Za, and Za/Z v indicates that the BLM prescription does indeed 
re-sum an important class of higher-order contributions. On the other hand, the coupling 
a Ms((/Ms) seems ' empirically, to work less well. In continuum perturbative QCD, it usually 
does not matter whether one adopts cty^qy), «ms(^ms") or some other renormalized cou- 
pling (at the BLM scale), once two- loop effects are included. It would not be surprising for 
the same to hold for short- distance quantities in lattice gauge theory, such as improvement 
coefficients. 
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VI. CONCLUSIONS 



In this paper we have compared non-perturbative calculations of several improvement 
coefficients to perturbation theory with the BLM prescription. Previously this could not be 
done, because the "BLM numerators" in Eqs. (27)-(p3|) were not available. We find that, 
for the scale-independent quantities considered here, the integration of the log fc 2 -weighted 
integrals is numerically straightforward. 

BLM perturbation theory for the current normalization factors Zj agrees very well with 
non-perturbative calculations of the same quantities. Here the leading power correction is 
only of order (Aa) 2 , and the small deviations can probably be removed with a two-loop cal- 
culation. Note that generalizations of the BLM method for higher-order perturbation theory 
have been considered in continuum perturbative QCD |2B| and in lattice gauge theory |37| . 

For the improvement coefficients bj and cj, the leading power corrections are of order 
Aa (and in the Schrodinger functional also of order a/L = 1/8), while some of the one-loop 
coefficients are small. It is consequently difficult to diagnose the discrepancies. By noting the 
size and dependence on g 2 of the differences, we concur with the authors of Refs. |24|, [25] , 
namely, that power corrections contaminate the non-perturbative results. In particular, 
it seems unlikely that higher orders in perturbative series could explain all discrepancies 
between one-loop BLM perturbation theory and the results from Refs. plL p2|. E3] . 
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